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In connection with the experiments recently achieved on doped crystals, biological samples, doped 
optical fibers and semiconductor heterostructures, we revisit the theory of the propagation of a 



pulse-modulated light in a saturable absorber. Explicit analytical expressions of the transmitted 
pulse are obtained, enabling us to determine the parameters optimizing the time-delay of the trans- 
mitted pulse with respect to the incident pulse. We finally compare the maximum fractional delay 
or figure of merit so attainable to those which have been actually demonstrated in the experiments. 
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I. INTRODUCTION 



Dynamics of saturable absorbers is often well reproduced by using a two-level model 
with a coherence relaxation-time very short compared to the population relaxation-time. 
The propagation of laser pulses in the medium is then simply described by two equations 
coupling the light intensity and the difference of populations. As far back as 1965, Gires 
and Combaud [lj used this model to analyze the transmission of laser pulses through dye 
solutions. They considered pulses of duration long compared to the population relaxation 
time, but this approximation is relaxed in subsequent works . Calculations made 

by Selden in this more general case enabled him to explain not only the narrowing of the 
transmitted pulse but also its skewing and the time-delay of its maximum [5|. Selden also 
studied the transmission of a laser beam when its intensity is slightly modulated by a low 



frequency sine- wave [6fl. He showed that the effect of the saturable absorber is to increase 
the modulation depth and to introduce aphase delay of this modulation. The experimental 
data obtained by Hillman et al. on ruby [7j are in full agreement with his predictions on the 
modulation depth. Although often overlooked, the above mentioned theories 
are applicable to most of the recent experiments achieved on various saturable absorbers, 



including doped crystals [8, 



wells 
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17| . quantum dots 



0, 
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13fl, biological film and solution 
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jl quant urn 
22|- Developed 



201 1 and doped optical fibers 
to attain pulse velocities as slow as possible, these experiments are currently analyzed in 
terms of coherent population oscillations (CPO), homogeneous hole-burning [23|| and group 



velocity. As extensively discussed in 



24 



261 ] . such an analysis is questionable. In most 



cases [27|], the population oscillations are not created in the medium under the combined 



action of two independent coherent beams 



23| but results from the intensity modulation 



of a single incident beam. The phenomenon is thus insensitive to phase and frequency 
fluctuations of the optical field. The group velocity, attached to a given optical frequency 
and obtained by expressing that the phase of the field is stationary near this frequency, 
looses then its relevance. We also remark that the identification of the group velocity to 
the ratio of the medium thickness over the time-delay of the pulse maximum, often made 
in the literature, is incorrect. As a matter of fact, the saturable absorption and the CPO 
approaches are based on the same approximations, namely that the coherence relaxation 
time is infinitely short compared to the population relaxation time, the Rabi period and 
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the inverse of the deviations of the laser frequency from the line frequency. In the CPO 
approach, analytical results have only been obtained in the particular case of a weak sine- 



wave modulation, 
same results 



13, 



25 



he fact that the saturable absorption approach then gives exactly the 



28| shows that the two approaches are equivalent. However the saturable 
absorption approach is more straightforward (it avoids the passage through the refractive 
index) and, as shown in the following, is more efficient since it provides analytical results in 
much more general situations, in particular not only when the pulse acts as a probe but also 
when its interaction with the medium is fully nonlinear. Finally the saturable absorption 
approach better corresponds to the experimental conditions where the inverse of the pulse 
duration is generally much smaller than the fluctuations of the optical carrier frequency. 

For the first time to our knowledge, we provide in the present paper explicit analytical 
expressions of the transmitted pulse with a special attention paid to its delay with respect to 
the incident pulse and to the optimization of this delay. In Section HD, we recall the general 
equations describing the propagation of intensity-modulated light in a saturable absorber. 
The case of pulses superimposed to a continuous background with a small modulation index 
is examined in Section HTT1 The nonlinear propagation of pulses in the absence of background 
and the general case (pulses and background of arbitrary intensity) are respectively studied 
in Sections HVl and IVl We finally compare in Section [VI] the fractional delays attainable with 
saturable absorbers to those which have been actually demonstrated. 

II. GENERAL ANALYSIS 

We consider a resonant light beam propagating in the z-direction through a saturable 
absorber modeled as a two level system. As indicated before, we assume that the coherence 
relaxation time is infinitely short compared to the population relaxation time, the Rabi 
period and the inverse of the deviations of the laser frequency from the line frequency. It is 



then possible to adiabatically eliminate thejpo 
order to obtain the two coupled equations 



arization in the Bloch-Maxwell equations in 



;_poiarizatic 



9N 

r— = N(l + I)-l (1) 
I - -aIN (2) 
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In these expressions, r is the population relaxation time, N is the population difference 
normalized to its value at equilibrium, t is the time retarded by the propagation time in the 
host medium (negligible compared to the delays considered in the following), / is the beam 



intensity normalized to the saturation intensity 



291 ] and a is the absorption coefficient in 



^(4 •/•!"/• '> --;• =0 (3) 



the linear regime. Combining Eqffl and Eqj2] we easily get the nonlinear wave equation: 

d_f m 

dz Ydt 
and the transmission equation 

T ~Et~ + Iout + ln Iout + aL = T ~§f + Iin + ln Iin ( 4 ) 



where L is the absorber thickness and I ou t is the normalized intensity of the output 

verys 

I out + ln lout + otL = I in + ln I in (5) 



(input) wave. When the input intensity is constant or verv slowly varying at the scale of r, 
EqfJ] is reduced to the well-known saturation equation 

mm 



Although established with a two level model, this equation fits very well the transmission 
curve of multilevel saturable absorbers. This result is illustrated FigJT] where we compare 



the predicted transmission to that actually measured on a erbium-doped optical fiber 



30]. 



III. CASE OF SMALL MODULATION INDEX 

We consider first the important case where the pulses (containing the useful signal) are 
superimposed to a large dc background C. The input and output intensities respectively 
read I in (t) = C in + s in (t) with s in (t) < C in and I out (t) = C ou t + s ou t(t) with s ou t(t) < C ou t- 
Making a calculation at the first order in s in {t) and s out (t) and taking into account Eqj5] 
relating C out and Cj n , we get : 

ds ut , Sout C ut ( dSin Si n \ . . 

+ — = ~Fi \—JT + — ( 6 ) 



dt T b C in \ dt T a 

where r a = 1/ (1 + Cj n ) and u = 1/ (1 + C ou t)- Assuming that Sj„(— oo) = 0, the general 
solution of Eq]6] may be written: 

t 



Sout(t) 



C ut 
Cin 



Sin{t) + (y-^rj e ' t/n / Sin{d> 9/n d6 



(7) 
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Figure 1: Transmission at 1530 nm of a erbium-doped optical fiber as a function of the incident 
power. The parameters are L = 7.5 m and aL = 6.8 . The points are experimental [3fJ] and the 
continuous line is obtained from Eq. [5] by adjusting the saturation power P sa t- The best fit is 
obtained for P sa t = —7.30 dBm , that is P sa t = 0.186 mW. The erbium concentration is small 
enough in order that energy transfer upconversion is negligible and that the absorption is fully 
saturable. 

I — I 

The impulse response h{t) [3j| is obtained by taking Si n (t) = 5(t) where 5(t) is the Dirac 
function. We get: 



h(t) 



a 



out 



Ci n 



5(t) 



n 



(8) 



where U(t) is the unit step function. Finally the transfer function [3jJ, Fourier transform of 
h(t), reads: 

H(Q) = ^(l+ _ ,?7. r A.^ = C - M ^ + l ^l (9) 



Ci, 



C in r a (1 + iVtr b ) 



The latter result can also be directly derived from Eqj6] by taking Si n {t) cc e int [25|| and 
is obviously applicable to the particular case of a sine-wave modulation, often used in the 
experiments. It is consistent with the previous calculations made in this case 0, [3| and 
with the experimental results. The phase delay of the intensity-modulation introduced by 
the medium has a maximum A$ m = tan -1 [(r a — u) /2^/r a Tb\ for Q = l/y/r a Tb. We remark 
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that A$ m < 7r/2, the upper limit being approached for Ci n 3> 1 and C out <C 1. Consequently 
the time-delay td of the output modulation can never exceed 25% of the modulation period 
T. 

Strictly speaking a sine- wave does not contain any information and, e.g., the previous 
delay td may also be seen as an advance T — td. An unambiguous demonstration of delay 
(or advance) requires to use pulses of finite duration and energy. Ultraslow "velocities" L/td 
can be achieved by using dense media with long relaxation times Q, [l^|. However, in view 
of potential applications, the important issue is not merely to achieve ultraslow light but to 
produce delays as large as possible compared to the duration of both the input and the output 
pulses. In the following we will thus characterize the slow light systems by their figure of 
merit or generalized fractional delay 

F = t d / max(r in , r out ) (10) 

where Tj n (T out ) is the full width at half maximum of the input (output) pulse. Our definition 
is identical to the usual one (F = t d /r in ) when r out « r in or r out < r in . 

We consider input pulses such that Sj n (i) is continuous, bell-shaped, symmetric and 
centered at t = 0. General properties of the output pulse s out (t) can be derived from the 
relations s out (t) = h{t) ® s in {t) or S out (Q) = H(Q)S in (Q) where S out (Q) and S in (Q) are 
respectively the Fourier transforms of s out (t) and s in (t). Since s in (t) is centered at t = 0, 
the center-of-mass t cm of s ou t(t) coincides with that of h(t). We get 

f +oc th(t)dt 

t = ^22 li = Ti-T 

Cm J^h(t)dt U ^ [ ' 

Since r a — rj (1 + Ci n ) and u = rj (1 + C out ), t cm will be always smaller than the relaxation 
time r, this limit being only attained when C in ^> 1 and C out 1. The largest pulse- 
delays are expected in these conditions but one should remark that, due to the distortion 
(asymmetric broadening), the delay td of the pulse maximum may strongly differ from t cm . 
Eqs. EM] show that the distortion will be negligible when r in ^> u (long pulses). We have 
then td ~ t cm and thus F <^ 1. More generally td will be as large as possible if the first term 
of EqJT] (not delayed) is small compared to the second one, that is again when C- m 3> 1 , 
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C out < 1, and thus t cm m r. We then get s out {t) m C out g(t) with 



g(t) 



T 



<8> s in (£) 



-t/T 



t (0)e e/T ^ 



FT' 1 



1 + iVLt 



(12) 



where FT~ l is a shorthand notation of the inverse Fourier transform. When Xj n <C r (short 
pulses), EqJT2l shows that td = 0(r in ) while r out ~ rln2 (the duration of £/(t)e~'/ r ) and 
thus F <C 1, as previously. A maximum of the fractional delay is expected for r in = 0(r) 
but its determination obviously requires to specify the pulse shape. 

We consider first the realistic case of pulses having a strictly finite duration (hereafter 
cos-pulses), such that s, n (t) = A in cos 2 (nt/2Ti n ) for — r in < t < r in and Si n (t) = elsewhere 
(FigEJ). Eqrjjthen leads to 



9(t) 



A: 



COS 



1 + 



M. i _|_ JUL s i n i jlL i _ | jzei. i P -(t+nn)/ 



(13) 



(W r ») 2 + 1 

for -T in <t< r in , g(t < -T in ) = and g(t > T in ) = g(r i „)e~ ( *~ T " l)/T '. As expected, s out {t) 
has an exponential fall at the end of the input pulse (t > T in ). The time-delay td of the 
maximum is given by the implicit equation: 



sin 



cos 



+ e 



— (*<J+7"in)/T 



(14) 



Asymptotic calculations show that td ~ t (1 — 7r 2 r 2 /2r| l ) for Tj n ^> r and that ~ 
Tj„ ^1 — 2t^ 2 /ht 1 / 2 ^ for r in r . When T/T in varies from to oo, t d /T in increases from 
to 1 while T out increases from T in to oo {r out rin 2 , see above). Starting from 0, the 
fractional delay F, equal here to t d /T out , begins to increase before to decrease to 0, in agree- 
ment with our general predictions (see inset of Figj2|). It attains its maximum F max = 31% 
for t/t^ = 0.9. This maximum is very flat since F max > 29% for 0.6 < r/rj n < 1.5. Figj2] 
shows the intensity profiles of the output pulses obtained for r/r in = 0.2, 0.9 and 5. 

Similar results are obtained in the classical case of gaussian pulses. Taking s in (t) = 



A in exp (-£ 2 / r p) with t p = r i „/2v / m2, we get: 



g(t) « Ai 



Tpy/TC 
" 2T 



1 + erf 



U 2t) 



t 



- 2 \ 1 



(15) 



where erf(x) is the error function. The optimal r/r^ (1.05) is close to that obtained with 
cos-pulses and F max is nearly the same (29%). The main difference is that the delay td is 
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Figure 2: Intensity profile of the output pulses obtained in the case of small modulation index for 
T/Vin = (a) 0.2 (b) 0.9 and (c) 5. The profile of the input pulse (cos-pulse) is given for reference 
(dashed line). The time unit is its full width at half maximum Tj n . Inset: Fractional delay as a 
function of the ratio r/rj n . 

no longer limited by r in . Delays td > r in can be obtained when r/r in 3> 1. Asymptotic 

In 



calculations then shows that td = r„ 



A delay td ~ r in is attained for 



T/r in ps 17. The output pulse is then very broad {r out ~ 12rj n and F ps 8% ). When 
the double condition Cj n ^> 1 and C out «C 1 is not met, the term proportional to Si n (t) in 
s out{t) ( see EqJT]) is not negligible and r b < r. The fractional delay is reduced accordingly. 
Considering, e.g., cos-pulses with C in = 1 and C out = 1/10 (attained by taking aL ps 3.2), 
we find F max ps 9% instead of 31% in the ideal case. 
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IV. PULSES WITHOUT BACKGROUND 



We consider now the case where C« n = 0, without restriction on the pulse amplitude. The 
medium being initially at equilibrium (N(— oo) = 1 ), Eqffl and Eqj2] show that N(t) > 
and s out (t) < s in (t) at every time. If the input pulse has a strictly finite duration (as the 
cos-pulses), s ou t(t) will thus stop at the same time that Si n (t). This result strongly contrasts 
with that obtained in the previous section (see Figj2]). 

When the input pulse is very short {r in <C r ), the population difference cannot follow 
the rapid change of the intensity and, roughly speaking, retains its initial value (sudden 
approximation). From Eqj2], we then retrieve the result corresponding to the linear regime, 
namely s out {t) = exp (—aL) s in (t). The pulse is only attenuated (neither distorted nor 
delayed). Conversely, when the input pulse is very long, s out {t) and Si n {t) are related by 
EqJHl The output pulse remains symmetric and centered at t = (no delay) but may be 
strongly narrowed |5fl. Finally, when Tj n and r are comparable, the output pulse will be at 
once narrowed, delayed and skewed. To study the general case, we consider the function 
Z(t) introduced by Selden Q| : 

Z(t) = In s out (t) - In s in (t) + aL (16) 
The transmission equation (EqHJ) then reads: 

+ Z = s m (t) [1 + e (z " aL) ] = s m (t) - s out (t) (17) 

with the initial condition Z{— oo) = 0. For given Sj n (t), EqJT7l shows that Z(t) and thus 
the shape of the output pulse will be independent of the optical thickness aL as early as 
the latter is large enough in order that s out {t) 1 and s out (t) s in (t) at every time. The 
pulse delay is expected to have then attained its maximum. We have checked this point 
by numerically solving EqJT7l Since we are mainly interested in maximizing the fractional 
delay, we will assume in the following that the previous condition on aL is actually met. 
EqJT7l is then reduced to: 

dZ 

r— + Z = s in (t) (18) 
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with the analytical solutions 



Z(t) 



e 



T 



I 



s m {6)e- e/T dO 



(19) 



— oo 



Soutit) = e 




(20) 



We see that Z(t) = g(t), where g{t) is the function introduced in Sec lIIII (Eqs fT2l [TBI and 
fTSi l. Consequently the delays t d considered in SeclIIII are now the delays tz of the maximum 
of Z{t) and thus of s out (t)/si n (t). Moreover, Si n (t) being centered at t — 0, Eqj20l shows that 
the new delay td of the pulse maximum will be smaller than tz and that, for an input pulse 
of given shape, t d (r out ) will be the larger (smaller), the larger is the amplitude A in . 

For a given amplitude A in , the shape of the output pulse and the fractional delay 
F = td/Ti n only depends on the ratio r/r in . For long pulses (r in 3> r ), EqJTBl takes the 
approximate form Z(t + r) ~ s in (t). We then get t z ~ r, s out (t) « e~ aL s in (t) exp [s in {t — r)) 
and, since ds in /dt = for t = 0, t^/r « + 1). For r/r in —* 0, F — ► as expected 

and s out (t) tends to the value given by Eqj5] so long as aL is actually large enough in or- 
der that s ou t{t) <^ 1 . Conversely when r/r in — > oo, dZ/dt — > 0, Z(t) — > Z(—oo) = , 
s ut{t) — > e~ aL Sj n (t) (as in the general case) and, again, F — > 0. Finally, a maximum of F 
(increasing function of A in ) will be obtained for an intermediate value of r/r in . 

Figure [3] shows the intensity-profiles of the output pulse obtained with cos-pulses for 
A in = 1 , 10 and 100 (keep in mind that A in is the peak intensity of the input pulse 
normalized to the saturation intensity). For each A in , r/T in is optimized in order to lead 
F to its maximum F max . Note that the narrowing of the output pulses is significant but 
that their skewing is moderate (fall steeper than the rise). We have systematically explored 
how F max , the corresponding T out /T in and tz/T in depend on the saturation for A in ranging 
from 0.2 to 10000 (FigfJ)). Since s in (t) stops at t = r in , the fractional delay cannot exceed 
unity. In fact, the limit F max = 1 is very slowly approached for very large values of A in . 
Asymptotic calculations then show that F max ~ 1 — (128 /-K A A in ) 1 ^ , this maximum being 
attained for r/r in w (2Af n /7r 2 ) 1//5 . Even for A in as large as 10000, F max is only 0.83. 

Comparable results are obtained with gaussian pulses for reasonable peak intensities of 
the input pulse, say for A in < 50 (Figj5]). For larger A in , some differences appear because 
the gaussian pulses have infinite wings. There is thus no theoretical limit to tz and t d . For 
example, F max slightly larger than 1 is attained for A in = 10000. 



10 




Figure 3: Intensity profile of the output pulses obtained in the case of an input pulse without 
background for Ai n = (a) 1, (b) 10 and (c) 100 with r/rj n = (a) 0.6 (b) 1.5 and (c) 4.2 (the value 
maximizing the fractional delay in each case). The profile of the input pulse (cos-pulse) is given for 
reference (dashed line). 

At this point, one should recall that the previous fractional delays F max will be actually 
attained only if the optical thickness is large enough in order that s out (t) <C s in (t) at every 
time, that is if exp [Z max — aL] <C 1 where Z max = Z{t z ). This condition is satisfactorily 
met for aL = Z max + 3. When A in is small (large), the optimum T/r in is also small (large). 
In the first case Z(t + r) w Sj„(t) (see above) and Z max ps A in C 1. In the second one, 
we easily get the asymptotic forms Z max ps rA in T in /T with r = 1 for cos-pulses and r = 
y / 7r/2v / ln2 ps 1.06 for gaussian pulses (1 Z max <^ A in ). For intermediate values of A in , 
Zmax < mm (An> r ^in T in/ T ) an d the condition s out (t) <^ s in (t) will be met in every case by 
taking 

aL = mm(A min , rA min r in / r) + 3 (21) 
Provided that r/T in is actually optimized to attain F max , the second condition of validity of 
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Figure 4: F max , tz/Ti n and T ou t/ri n as functions of the peak intensity Ai n in the case of cos-pulses. 
For Ai n > 1000 , the ratio r/Tj n maximizing F and F max itself are well approximated by the 
asymptotic formula r/T in « {2A 2 in /-K 2 ) 1 ^ and F maa . « 1 - (128/7r 4 J 4 in ) 1//5 . 

our calculation, namely s out {t) <C 1, is then automatically fulfilled. 

V. PULSE AND BACKGROUND OF ARBITRARY INTENSITY 

Comparing the results obtained with input pulses superimposed to a large background 
(Sec JIIII) and with pulses without background (Sec lIVI) . we see that the former are broadened 
in the medium with a rise significantly steeper than the fall (FigGJ) whereas the latter are 
narrowed with a fall steeper than the rise (Figj3|). We may then hope that better results 
will be obtained by using pulses superimposed to a suitably adjusted background. We 
thus consider in this section the case where Unit) = C in + s in (t) without restriction on the 
amplitudes of Ci n and Sj n (i). As previously and for the same reasons, we assume that aL 
is large enough in order that I ut(t) 1 and I ut(t) hn{t) at every time. By redefining 
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Figure 5: Same as Fig. [4] in the case of gaussian pulses. The ratio t/t^ maximizing F is 0.6, 1.5, 
5.0, 14 and 42 respectively for A in = 1, 10, 100, 1000 and 10000. 

Z(t) as Z(t) = \iaI out {t) — \ial in (t) + aL — Cj n , we find that EqH8] is unchanged and thus 
that Z(t) = g(t) as previously. In other respects the new definition of Z(t) leads to 

C out + s out (t) = [C m + s in (t)\ e^W-^ (22) 

Since s in (t), s out (t) and g{t) cancel for t = ±oo, C out = C in exp [C in — aL) in agreement 
with Eqj5]in the limit C out <C 1 considered here. Finally s out {t) reads 

Soutit) = [C m (e*<*> - 1) + s m (t)e 9it) ] e^" +Qi ) (23) 

When Ci n = 0, we retrieve the result given in the previous section (EqEZOl). Conversely when 
the modulation index is small, e 9 ^ — 1 ~ g{t and we get 

Soutit) = [s m {t) + C m g{t)}^ (24) 

'—'in 

a result consistent with EqJTl again in the limit C out <C 1 where r b « r. 
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Eq.23 enables us to determine the profiles of the output pulses for arbitrary values of 
the ratio Ci n /A in . We give FigJB] different profiles obtained when the input peak-intensity 
is fixed (Ci n + A in = 10). For each value of Ci n /A in , r/T in has been optimized in order to 
maximize F. As expected the addition of a background widens the output pulse. It does 
not significantly enlarge the attainable fractional delay, which very slightly increases as a 
function of Ci n /A in before falling down to the value calculated in the small modulation-index 
limit (see inset of Figj6)). However we remark that the resemblance of the output pulse to the 
input one can be improved by the presence of a background (see the profile of the output 
pulse obtained for Ci n /A n =0.54 . The latter effect has been recently demonstrated in a 



saturable gain system 



321 ] . The qualitative behavior shown Figj6] is general and is observed 



for any bell-shaped input pulse. 



VI. SUMMARY AND DISCUSSION 

We have theoretically studied the transmission of a pulse-modulated light in a saturable 
medium modeled as an ensemble of 2-level atoms with a coherence relaxation time extremely 
short compared to the population relaxation time. This model of saturable absorber gives 
theoretical results in good agreement with the experimental results obtained in the currently 
called CPO based slow light experiments. This was already pointed out in Ref.[3| about 
the representative experiments achieved on ruby Q|, Er 3+ :Y 2 Si05 crystal biological 
bacteriorhodopsin 



U] and quantum dots [18||. We checked that it is also true for the 



extensive experiments recently realized on erbium-doped optical fibers [22j. More specifically, 
we verified that, except for ultrahighly doped fibers (ion density exceeding 3 x 10 25 m~ 3 ), 
the maximum phase delays A$ m attained for a sine- wave modulation and the corresponding 
modulation frequency are in agreement with those given by the model (see the discussion 
following Eqj9]). 

Thanks to the relative simplicity of the transmission equation of the model system (Eqfjj), 
it has been possible to obtain explicit analytical expressions of the output pulse and to 
optimize the figure of merit or fractional delay F of the system (Eq JTOj) . Our main findings 
are as follows. When the input pulse sits on a much larger dc background Cj n (intensity 
normalized to the saturation intensity), the output pulse is asymmetrically widened with a 
rise steeper than the fall. This behavior is qualitatively analogue to that of the usual slow 
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Fractional Time t/x in 

Figure 6: Intensity profile of the output pulses obtained for an overall peak intensity Ai n + Ci n = 10. 
In order of increasing width, the represented profiles correspond to Ci n /Ai n = 0, 0.11, 0.54 and 
9.0. For each value of Ci n /Ai n , r/Ti n is optimized in order to maximize the fractional delay. The 
profile of the input pulse (cos-pulse) is given for reference (dashed line). Inset: F max as a function 
of C in /A in ■ 

light systems (see, e.g., j^]) but the pulse shape may be much more asymmetric, with an 
exponential or nearly exponential fall (Figj2]). The fractional delay F depends on Cj n , on 
the linear optical thickness aL (which determines the intensity of the output dc background 
C out ) and on the ratio T/r in of the population relaxation time over the width of the incident 
pulse. It attains its maximum F max ps 30% (slightly depending on the precise shape of the 
input pulse) when C in ^> 1, C out <C 1 and r/r in « 1. When C in = 1 and C out = 1/10 
(aL « 3.2), F max falls down to 9%. Larger fractional delays are obtained by using input 
pulses of large peak intensity A in without background. Contrary to the previous case the 
output pulse is now narrowed with a fall moderately steeper than the rise (FiglBj). The 
largest fractional delays are attained when A in is as large as possible (Figs. I4f5l) provided 
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that the optical thickness is itself very large (Eqj2TT). Note that the ratio r/r in maximizing 
F also increases with A in . In the reference case A in = 10, F max « 36% for r/r in « 1.5 and 
aL « 10.7 [3jj|. Finally, for a fixed value of the overall peak intensity of the input beam, 
the addition of a dc background does not significantly enhance the fractional delay but may 
improve the symmetry of the output pulse (Figj6]). 



In fact, there are few time-reso 
evidence of pulse delays [8 
Eq JTOl) are all smaller than 



10 



ved experiments on saturable absorbers giving direct 



1&I21 



|22j. The obtained fractional delays (as defined 
o. There are different reasons for that. The main one is that 
the input intensities Ci n and/or A in are too small, typically of the order of 1, at the best of a 
few units. Second the linear optical thickness is not adapted. Third the pulse duration is not 
optimized. The erbium-doped optical fiber seems a good candidate for the demonstration of 
a larger fractional delay. The saturation power is low (< 0.5mW) and normalized intensities 
C in and/or A in of 100 can be easily achieved. A fractional delay of about 60% (Figj3b) 
would then be attained with an input pulse of duration r in rs 0.23r rs 2.4 ms and a linear 
optical thickness aL « 23. The latter would be obtained in a fiber of reasonable length 
(L < 4 m) with an ion density p ~ 2 x 10 25 m -3 22J for which the saturation model is valid. 
Note that larger fractional delays (up to 1.5 with our definition) have been demonstrated in 
undoped fibers by exploiting Brillouin scattering (35} but this result is obtained with much 
longer fibers. 

We finally remark that the pulse-delay mechanisms in a saturable absorber strongly differ 



from those involved in the "pure" slow-light experiments |36| . The former are nonlinear 
and non coherent whereas the latter are linear and coherent. Moreover the propagation 
phenomena are essential in the second case whereas they are absent in the first one. This 
point is illustrated by our calculations made for an input pulse of strictly finite duration 
(Sec. HVl) . We have shown that the output pulse then stops at the same time that the input 



one. On the contrary the propagation effects are responsible of an important delay in the 



linear case. This explains in particu 
with an electromagnetically induced 



ar the very large fractional delays attained in media 



37|] or a natural [38|, |39|] transparency window. 
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